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Percolation model for structural phase transitions in Lii_ x H x I03 mixed crystals 

P. H. L. Martins 1 ' 2 *, J. A. Plascak 2 t, and M. A. Pimenta 2 * 
1 Universidade Federal do Tocantins, Caixa Postal 111, 77001-970, Palmas, TO - Brazil 
2 Departamento de Fisica, Instituto de Ciencias Exatas, Universidade Federal de Minas Gerais, 
Caixa Postal 702, 30123-970, Belo Horizonte, MG - Brazil 
(Dated: February 2, 2008) 

A percolation model is proposed to explain the structural phase transitions found in Lii-^H^IOs 
mixed crystals as a function of the concentration parameter x. The percolation thresholds are 
obtained from Monte Carlo simulations on the specific lattices occupied by lithium atoms and 
hydrogen bonds. The theoretical results strongly suggest that percolating lithium vacancies and 
hydrogen bonds are indeed responsible for the solid solution observed in the experimental range 
0.22 < x < 0.36. 
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The lithium iodate (a-LiI03, hexagonal lattice) and 
the iodic acid (a-HI03, orthorrombic lattice) have been 
both intensively investigated, owing to their interest- 
ing properties and useful optical applications^ 2 * 3 ^. The 
lithium iodate is also interesting for fundamental rea- 
sons, specially those related to its structural phase 
transitions 5 . It is known that most of the properties of 
the a-LiI03 are strongly affected by its growth condi- 
tions, related to three mainly factors: the temperature, 
the evaporation rate, and the pH of the mother LHO3- 
H2O-HIO3 solution, which controls the quantity of im- 
purities in the lattice. In this case, the impurities can 
be hydrogen atoms. Mixed crystals are easier grown us- 
ing more acid solutions (with a greater HIO3 concentra- 
tion), but they have worse optical quality. hii- x H x l03 
mixed crystals are also important for technical applica- 
tions, since they present a high piezoelectric coefficient 6 . 

Although the structures of iodic acid and lithium io- 
date are known for more than sixty years, the structure 
of the mixed lithium iodate-iodic acid solid solution has 
been refined only in the last few years**"^. Regarding 
the lithium substitution, neutron diffraction studies have 
shown that the hydrogen does not substitute the posi- 
tion of the lithium. Instead, when lithium goes out of 
the lattice, its site becomes empty, while the hydrogen 
forms hydrogen bonds between oxygens that belong to 
different iodate groups^. 

Ricci and Amronifi first mentioned the existence of the 
so-called solid solution of Lii-^H^IOs for x varying con- 
tinuously over the range 0.22 to 0.36. Since then, much 
effort has been dedicated in studying these compounds, 
but the reason of these particular values is not yet well 
understood. In this work we propose a percolation model 
that explains the lower and upper limits experimentally 
observed for crystallizing the solid solution of these mixed 
crystals. 
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Let us first discuss the principal properties of pure 
lithium iodate. Structural phase transitions on this crys- 
tal have been largely studied 5-11 - 12 - 13 . At room tempera- 
ture and low pressure, two forms are possible: a (hexago- 
nal) and f3 (tetragonal) . The a structure has been known 
since 1931— and is the more stable one at room tem- 
perature. There are two molecules per unit cell with 
the lithium atoms disposed on a simple hexagonal lat- 
tice. Lattice constants in A, with errors in parenthe- 
sis, are given by 3 : a = 5.48(1) and c = 5.18(1). The 
Li-0 distance is approximately 2.12 A. As the temper- 
ature increases, a-LiI03 exhibits two first-order phase 
transitions: a — > 7 — > j3, where phases 7 and ft have 
orthorhombic and tetragonal symmetries, respectively^. 
Several works reveal discrepancies in the transition tem- 
peratures as well as in the temperature range where the 7 
phase exists. These differences are attributed to, among 
other reasons, the quantity of HIO3 impurities in the 
lattice**^ 2 *— Thus, a good understanding of the al- 
terations caused by the presence of hydrogen atoms is 
essential for studying these phase transitions. 

Let us now turn to the a-HI03 crystal. The structure 
of a-iodic acid, though not as simple as that of the a- 
LHO3, has been known since 1941, by means of X-ray 
diffraction—. It is orthorhombic, with four molecules per 
unit cell and lattice constants in A: a = 5,520(5); b — 
5,855(5) and c = 7,715(5). The iodate group is pyra- 
midal, although slightly distorted, and has the following 
interatomic distances in A(with error of 0.04 A or less): 
I Oi = 1.81; I-On = 1.89; I-O m = 1.80; Oi-O n = 2.75; 
Oi-O m = 2.78 and On-Om = 2.78, where Oi, O n and 
Oin are the three oxygens of the iodate group. The hy- 
drogen is bonded to the oxygen On, at a distance of 
1.01 A. This oxygen has two other neighbors: Oj and 
Ojjj. Here, prime and double primes denote different io- 
date groups. On forms two hydrogen bonds, On-H-OJ 
and On-H-Ojjj, of equal intensities. The distances be- 
tween the two oxygens indicate that the hydrogen bond is 
strong (there are here two bonds per iodate group each 
of which is about the strenght of the single hydrogen 
bond in water and ice)**V A better way to visualize this 
structure is to look at a projection onto the vz plane, 
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as shown schematicaly in Fig. ^ where the whole IO3 
group has been represented by just one symbol (and no 
distinction is made among the three oxygens). There ex- 
ists a kind of bifurcated hydrogen bond: a single proton 
is able to form two hydrogen bonds, i.e., the hydrogen 
closest to an iodate forms a bifurcated bond with two 
other iodates. Thus, there are two hydrogen bonds for 
each iodate group. These bifurcated bonds and the weak 
bonds between iodine and oxygens belonging to other io- 
date groups hold the iodates together. 
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FIG. 1: Schematic representation of the HIO3 structure onto 
the yz plane. Bigger circles represent IO3 groups and the 
smaller ones hydrogen atoms. Bifurcated hydrogen bonds are 
represented by solid lines which bifurcate on dashed lines at 
the hydrogen atoms (after Ref^-). 

The existence of bifurcated hydrogen bonds has been 
known since 1939, when the structure of glycine was 
determined^. Nevertheless, with the exception of a few 
early studiesiii2*i^, this subject has not received much 
attention. The term bifurcated is used to denote two dif- 
ferent kinds of bond. In a bifurcated acceptor, two dif- 
ferent acceptor oxygens share the same donated proton. 
In the case of a bifurcated donor a single oxygen donates 
two protons to another single oxygen. The first case is 
what occurs in the structure of iodic acid, as determined 
in Refil 

Several works have been devoted to the mixed crys- 
tal Lii^ x HJOj 7 i 10 i 20 i 21 i 22 i 23 . Ricci and Amronia first 
reported the existence of the solid solution for hydro- 
gen concentrations x in the range 0.22 < x < 0.36. In 
1990, Le Roy et ali concluded that, for x = 0.33, the 
mixed crystal possesses an average hexagonal symmetry, 
with an arrangement close to the HCP (hexagonal close- 
packed). The structure at room temperature is shown in 
Fig. [3 Indeed, they observed that protons do not sub- 
stitute the lithium atoms (2a in Wyckoff notation), but 



are randomly distributed through the lattice, occupying 
general positions (6c) between two oxygen atoms that 
belong to different iodates. Raman investigations^ are 
consistent with this structure and show that the random 
presence of hydrogens changes the selection rules for the 
Raman scattering for x in the range 0.27 < x < 0.36. The 
lattice constants are found to vary linearly with the con- 
centration x, according the relations'": a = 5. 464+0. 282a; 
and c = 5.165 — 0.622a;, with a and c in A. So, it is pos- 
sible to obtain x via: x — 5.985 — 6.335c/a, valid for 
0.22 < x < 0.36. 




FIG. 2: The crystal structure of Lii-zEUIOa (x ~ 0.33) pro- 
jected onto the xy plane, according Refi— . In this case, bifur- 
cated bonds are not considered. The protons do not replace 
the lithium ions, but are randomly distributed in general posi- 
tions between two oxygen atoms belonging to different iodate 
groups. 

As mentioned before, only in the range 0.22 < x < 0.36 
the solid solution is constituted as a single phase. Out 
of this range, there is a mixture of solid solution and the 
pure compounds (L1IO3 and HIO3). Thus, these limits 
seem to be favorable for crystallizing the mixed crystal. 
Our purpose is identifying these values as the percolation 
thresholds on the corresponding lattices. In this man- 
ner, we expect that 0.22 correspond to the critical point 
for site percolation on the LHO3 (simple hexagonal) lat- 
tice. Similarly, 0.36 might correspond to the percolation 
threshold for bond percolation on the HIO3 (orthorhom- 
bic) lattice. This seems to be plausible since, as has been 
discussed above, it has been observed that the protons 
do not substitute the lithium atoms, but are randomly 
distributed through the lattice, occupying general posi- 
tions between two oxygen atoms that belong to different 
iodate o 7,23 . The proton positions should obey the ar- 
rangement depicted in Fig. 1. As the proton mediates 
the hydrogen bond, the net effect is a random dilution of 
hydrogen bonds on the crystal structure of Fig. 1. Thus, 
this system turns out to be a quite interesting problem 
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involving site and bond dilution in the same crystalline 
structure. 

It is known that, when lithium atoms are removed from 
their sites on the original hexagonal lattice, vacancies are 
left on their positions^. As more lithium atoms go out 
(which corresponds to an increasing x) , clusters of neigh- 
boring vacancies are formed. There exists a given critical 
concentration x c at which the lithium vacancies percolate 
through the lattice. In other words, for vacancies concen- 
trations x such that x > x c , there are an infinite cluster 
of neighboring vacancies. This fact should means that, 
for x > x c , the presence of vacancies causes an instability 
on the LiI03 hexagonal lattice and then propitiates the 
crystallization of the solid solution. 

On the other hand, concomitantly to the removal of 
lithiums, hydrogen atoms are added in the lattice, orig- 
inating H-bonds. By an analogous argument, when the 
number of H-bonds is sufficiently large, there will be a 
cluster of neighboring bonds that percolates through the 
lattice, being possible to consolidate and stabilize the 
HIO3 orthorhombic lattice. 

We performed Monte Carlo simulations in order to ob- 
tain the percolation thresholds on the corresponding lat- 
tices. The algorithm we have used is that due to New- 
man and ZifS4, applied to three-dimensional lattices 2 ^. 
Despite to the fact that the percolation threshold for 
the hexagonal lattice need not be recomputed, since its 
value is quite well known from the literature [for instance, 
Pc = 0.2625(2)^, p c = 0.2623(2)2^, we have done just 
one simple simulation on a large lattice to see the per- 
formance of the method in order to apply it to the less 
studied case of the present bond dilution. In the site di- 
lution problem, for a given system size, we start with an 
empty lattice. So, we fill the lattice, choosing sites to 
be occupied, at random. When a new site is added, it 
can originate an isolated cluster (if all its neighbors are 
empty), or it joins together two or more clusters. After a 
complete filling of the lattice (which corresponds to one 
Monte Carlo step), we can evaluate quantities of inter- 
est, like the percolation threshold. Repeating this proce- 
dure several times, we get more accurate results. Using 
periodic boundary conditions, the percolation threshold 
is the concentration of occupied sites at which, for the 
first time, a cluster of neighboring occupied sites wraps 
around the system. The bond problem is treated in an 
analogous way. 

For simulating the removal of the Li atoms, we can 
think that it is necessary to begin with a full lattice and 
then take out sites at random. In this sense, percolation 
means that there exists an infinite cluster of vacancies. 
If we replace Li atoms by empty sites and vacancies by 
occupied sites, we restore the standard percolation prob- 
lem. 

We applied the cited algorithm on the simple hexago- 
nal lattices with up to 8 000 sites, which corresponds to 
a linear dimension L = 20. As we are interested only in 
the lithium positions, each site represents a lithium atom 
(the iodates localization are not relevant in this case). 



Performing 2.0 x 10 5 Monte Carlo steps, one obtains 
x c = 0.26 for the percolation threshold on this lattice. It 
is clear from the above result that for the purpose of the 
present study just this lattice gives a value which is quite 
comparable to those listed above from a more accurate 
finite-size scaling approach, namely, p c = 0.2625(2) 2 - 
and p c — 0.26 23(2) 26 . Within our conjecture, this value 
might correspond to the lower limit for crystallizing the 
solid solution of Li1_.rH.rIO3. 

Due to the complex structure of HIO3, some modifica- 
tions on the algorithm concerning the bond dilution are 
needed, in comparison with the usual percolation model. 
Since hydrogen bonds are formed between different io- 
date groups, we can represent, in our simulations, each 
iodate by a single site in the lattice. One has the restric- 
tion that bonds can exist only between sites that are lo- 
cated on different planes in order to reflect the fact that, 
in the real iodic acid, the distance between neighboring 
iodates on the same plane is greater than the maximum 
value for existing hydrogen bonds. The presence of bifur- 
cated bonds leads to the additional constraint that one 
has always to choose a couple of bonds together. Since 
we know a priori the bonds that share the same hydro- 
gen, we can choose a site (iodate) and then fill the two 
corresponding bonds. Using lattices with a total num- 
ber of sites up to 6 912 (which corresponds to a linear 
size L = 12, since there are four molecules per unit cell), 
and 5.0 x 10 5 Monte Carlo steps for each system size, 
simulations yield x c — 0.33 for the critical hydrogen con- 
centration on the specified crystal structure. Assuming 
the same accuracy as the previous site diluted model we 
believe that the number above is quite fine to be com- 
pared with the experimental one. A finite-size scaling 
estimate, providing a critical concentration with more 
digits of precision, would be irrelevant for the present 
purpose. 

Before achieving these results, other similar hypothe- 
ses were tested, although all of them sharing the same 
fundamental concepts. The first attempts did not in- 
clude bifurcated bonds. Bonds were choosen randomly 
one by one. Besides results were unsatisfactory, there is 
an incompatibility between the concentration of lithium 
vacancies and the corresponding concentration of hydro- 
gen bonds. This can be explained as follows. For x = 
(pure LHO3) there is one lithium for each iodate. Imag- 
ine a lattice with a total of N lithium atoms. So, there 
are also N iodates. Since each iodate has four neighbors 
there are 2N possible bonds. Removing one lithium, we 
need to introduce one hydrogen. If each hydrogen forms 
a single bond, when all N lithium atoms will be removed 
(x = 1), there will N hydrogen bonds, which corresponds 
to a fraction N/2N = 0.5 of the total number of bonds. 
This fact precludes the equivalence between the hydrogen 
concentration x and the related percolation thresholds. 
Bifurcated bonds avoid this problem, since N hydrogens 
form 2N bonds. Thus, the fraction of vacancies is equal 
to the fraction of bonds. 

According to our results, the solid solution should ex- 
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ist as a single phase for x in the range 0.26 < x < 0.33. 
In this sense, when x rises from to 0.26, lithium atoms 
are removed from the lattice but the a-LiI03 structure 
remains stable. For x > 0.26, the number of vacancies is 
so great that the a-LiI03 lattice does not support itself. 
Concomitantly, while x < 0.33, the hydrogen bonds are 
insufficient for supporting the HIO3 orthorhombic struc- 
ture. Thus, in the range 0.26 < x < 0.33 the system is 
disordered and we find the so-called solid solution. For 
x > 0.33 the percolating cluster of bonds maintains the 
HIO3 lattice stable. 

We conclude that our results are in good agreement 
with the experimental values and, although not stated 
in the experimental works, within the experimental er- 



rors. [For instance, the smaller concentration varies from 
0.23iS to 0.2'?2£]. Thus, our results evidence the perco- 
lating character on these mixed compounds and clearly 
indicate that the crystal stability is, in fact, related to the 
percolation thresholds supported by the corresponding 
crystal structures. Further experimental studies regard- 
ing the structural transition of this material would be 
very welcome, mainly concerning its second-order char- 
acter. 
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